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In this work we will focus on the effects produced by topological disorder 
on the electronic properties of a graphene plane. The presence of this type 
of disorder induces curvature in the samples of this material, making quite 
difficult the application of standard techniques of many body quantum theory. 
Once we understand the nature of these defects, we can apply ideas belonging 
to quantum field theory in curved space-time and extract information on 
physical properties that can be measured experimentally. 

1 Introduction 

Graphene is a two dimensional material formed by isolated layers of carbon 
atoms arranged in a honeycomb-like lattice. 

Each carbon atom is linked to three nearest neighbors due to the sp 2 hy- 
bridization process, which leads to three strong a bonds in a plane and a 
partially filled 7r bond, perpendicular to the plane. These 7r bonds will deter- 
mine the low energy electronic and transport properties of the system. 

It is possible to derive a long wavelength tight binding hamiltonian for the 
electrons in these 7r bonds([Wallace(1947)]). This hamiltonian is: 



where vf being a constant with dimensions of velocity (vf ~ 10 3 to/s). The 
wave equation derived from the hamiltonian (1) is the Dirac equation in two 
dimensions with the coefficients 7 J being an appropriate set of Dirac matrices. 
We can set for instance, 7 1 = 1 <g) u\ and -f 2 = T3 (g> 02, where the u, r matrices 
are related to the sublattice and Fermi point degrees of freedom respectively. 
The unexpected form of the tight-binding Hamiltonian comes from two special 
features of the honeycomb lattice: first, the unit cell contains two carbon atoms 
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belonging to different triangular sublattices, and second, in the neutral system 
at half filling, the Fermi surface reduces to two nonequivalent Fermi points. 
We will study the low energy states around any of these two Fermi points. The 
dispersion relation obtained from (1) is e(k) = ±?^|k|, leading to a constant 
density of states, p°(u>) = f 

2 A first model for the topological defects in graphene 

Several types of defects like vacancies, adatoms, complex boundaries, and 
structural or topological defects have been observed experimentally in the 
graphene lattice ([Hashimoto et al.(2004)Hashimoto, Suenaga, Gloter, Urita, and Iijima]) 

and studied theoretically (see for example [Lehtinen ct al.(2003)Lehtinen, Foster, Ayuela, Krasheninnikov, Nor 
[Lopcz-Sancho ct al.(2006)L6pcz-Sancho, Staubcr, Guinea, and Vozmediano]). 
Topological defects are produced by substitution of an hexagonal ring of the 
honeycomb lattice by an n-sided polygon with any n. Their presence impose 
non-trivial boundary conditions on the electron wave functions which are diffi- 
cult to handle. A proposal made in [Gonzalez et al.(1992)Gonzalez, Guinea, and Vozmediano] 
was to trade the boundary conditions imposed by pentagonal defects by the 
presence of appropriate gauge fields coupled to the spinor wave function. A 
generalization of this approach to include various topological defects was pre- 
sented in [Lammert and Crespi(2004)]. The strategy consists of determining 
the phase of the gauge field by parallel transporting the state in suitable form 
along a closed curve surrounding all the defects. 

&(6 = 0) = T c &(6 = 2tt) <s> V{8 = 0) = exp( I A a T a dv)y{6 = 2tt), (2) 

J c 

where A a are a set of gauge fields and T a a set of matrices related to the 
pseudospin degrees of freedom of the system. 

When dealing with multiple defects, we must consider a curve surrounding all 
of them, as the one sketched in figure (1): 

The contour C is made of closed circles enclosing each defect and straight 
paths linking all the contours to a fixed origin. The parallel transport operator 
Pc associated to the closed path is thus a composition of transport operators 
over each piece: 

P = P 7 i • Pi • P"! 1 • ... • P lN ■ P N ■ P" 1 . (3) 

As explained in [Lammert and Crespi(2004)]the total holonomy turns out 
to be 3 : 



3 The usual chiral lattice real vector basis for the honeycomb lattice is used in this 
derivation. 
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Fig. 1. Scheme of a prototypical curve enclosing multiple defects in which the state 
will be parallel transported. 

P=(i^ {T2 ) N exp(^(N+-N-)a 3 yxphfp^n j -m j )T)j . (4) 

From equation (4) we see that we have in principle three different gauge 
fields to incorporate into the Dirac equation, which couple to the matrices 03, 
T2, and T3 and whose associated fluxes are adjusted from (2). 



3 Generalization of the model 

In spite of its elegance, the model presented in the previous section does 

not contain the effects due to the curvature of the layer in the presence of 

these defects. The model can be generalized to account for curvature effects 

([Gonzalez et al.(1992)Gonzalez, Guinea, and Vozmediano], [Kochetov and Osipov(1999)]) 

by coupling the gauge theory obtained from the analysis of the holonomy in 

a curved space. 

The substitution of a hexagon by a polygon with n < 6 sides gives rise to 
a conical singularity with deficit angle (27r/6)(6 — n), which is similar to the 
singularity generated by a cosmic string in general relativity. The Dirac Equa- 
tion for a massless spinor in a curved spacetime is ([Birrell and Davis(1982)]): 

i^ix^-rf^^o, (5) 

(T) 

where J^- is a set of spin connections related to the pseudospin matrices in 
(4) and 7^(3;) are generalized Dirac matrices satisfying the anticommutation 
relations 
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{r(x),r(x)}=2g^(x). (6) 

The metric tensor in (6) corresponds to a curved spacetime generated by 
an arbitrary number of ./V parallel cosmic strings placed in (a,, b{) (here we will 
follow the formalism developed in [Aliev et al.(1997)Aliev, HOrtacsu, and Ozdemir]): 

ds 2 = -dt 2 + e- A ^' y \dx 2 + dy 2 ), (7) 

with A(x,y) = J2iLi ^Mi logdO^ — a i) 2 + ( x ~ hi) 2 ] 1 / 2 ). The parameters fa are 
related to the angle defect or surplus by the relationship a = 1 — Afa in such 
manner that if < 1(> 1) then fa > 0(< 0). 

From equation (5) we can write down the equation for the electron propagator, 
Sf(x, x 1 ): 

i^(x)(d^-r^)S F (x,x') = -L<5 3 (x-x'). (8) 

v-g 

The local density of states N(u), r) is obtained from the solution of (8) by 
Fourier transforming the time component and taking the limit r' — ► r: 

N(u, r) = ImTrS F (LU, r, r). (9) 

Provided that we only consider the presence of pentagons and heptagons, 
the parameters fa are all equal and small (fa = fi = 1/24). We will solve 
equation (8) perturbatively in /z. 

When dealing with equation (8) we will reduce the number of spin connections 
derived in the previous section by the following considerations: First, we will 
consider an scenario where the number of pentagonal and heptagonal defects 
is the same - so the total number of defects is even. This suppresses the con- 
tribution from the first exponential in (4) . If we consider that pentagonal and 
heptagonal defects come in pairs as usually happens in the observations, we 
can neglect the effect of mixing of the the two sublattices that each individual 
odd-sided ring produces and hence eliminate the spin connection related to r 2 
from (8). Furthermore, we can disregard the spin connection related to T3 by 
the following argument: We will solve equation (8) perturbatively to first or- 
der of the parameter /x. In general if S F is the unperturbed Dirac propagator 
and V(ui, r) the perturbation potential, the first term of such solutions is: 

Sfa, r, r') - ^ J d 2 v"S F {u, r, r")V>, t")S° f {u>, r", r'), (10) 

and we trace S F (u), r, r) in order to get the first contribution to the density of 
states SN(lu, r). The trace operation eliminates all the terms appearing in (10) 
which are proportional to a traceless matrix, including the matrix related to 
T2- In fact, up to this order in perturbation theory, the only term that survives 
will be the one proportional to 7 . With all this in mind, the relevant spin 
connection terms are: 
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Fig. 2. First order correction to the local density of states in a region around two 
pairs of heptagon-pentagon defects located out of the image for increasing values of 
the energy. 

A(r) - ~^d v A,r 2 {T) = —tV^A (11) 

After all these simplifications we can write equation (8) in a more suitable 
form. Expanding the terms in (11) in powers of /j, we get the potential V(lu, r): 

V(u, r) = -2A~/°uj + iAfdj + ^{djA). (12) 

As we said, expression (10) gives us the first correction to the local density of 
states in real space. In figure (2) we present an example of the results obtained. 
We show the first order correction to the local density of states coming from 
two pairs of heptagon-pentagon defects located out of the image for increasing 
values of the energy. What we see is that as the frequency increases, the local 
density of states is enhanced and inhomogeneous oscillations are observed in 
a wide area around the defects. The spatial extent of the correction is such 
that the relative intensity decays to ten percent in approximately 20 unit 
cells. The model described in this work can be applied to other configurations 
of defects, such as simple pairs or stone- Wales defects. These results can be 
found in ([Cortijo and Vozmediano(2006)]). 
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